Coherent manipulation of electron spins up to ambient temperatures in Cr 5+ (S=l/2) 

doped K 3 NbO s 
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We report coherent spin manipulation on Cr 5+ (S = 1/2, 1 = 0) doped KsNbOs, which constitutes 
a dilute two- level model relevant for use as a spin qubit. Rabi oscillations are observed for the first 
time in a spin system based on transition metal oxides up to room temperature. At liquid helium 
temperature the phase coherence relaxation time T2 reaches ~ 10 fjs and, with a Rabi frequency of 
20 MHz, yields a single qubit figure of merit Qm of about 500. This shows that a diluted ensemble 
of Cr 5+ (S = 1/2) doped KaNbOs is a potential candidate for solid-state quantum information 
processing. 



Recently, electron spins in solids have been intensively 
discussed in terms of physical implementations of a quan- 
tum computer [l[ and several proposals for embodying 
solid-state spin qubits have been put forward. The dis- 
cussed physical systems comprise quantum dots [1, H, Hj] , 
phosphorous donors in silicon [f|, endohedral fullerenes 
j(| , nitrogen- vacancy centers in diamond @, H, H. HI EH , 
molecular magnets [12, EH and rare-earth ions 14 1 . They 
commonly make use of the well-characterized discrete en- 
ergy levels arising from the spin, orbital, or charge states. 
Even though transition metal ions have these states, to 
our knowledge, they have not been exploited as a basic 
building block of solid-state spin qubits. More impor- 
tantly they could be made essentially free of magnetic 
anisotropy and therefore suitable for on-chip deposition 
and spin manipulation (vide infra). 

In this Letter, we explore this possibility by lightly 
doping a S = 1/2 Cr 5+ ion into the nonmagnetic matrix 
of H^NbOs- Cr was chosen since its dominant isotope 

Cr (90.5% natural abundance) has nuclear spin I = 0, 
thus obviating complications like spin decohcrence due to 
hyperfine interactions. A further advantage of this sys- 
tem is that, in principle, it can be isotopically enriched 
with 53 Cr (J = 3/2, 9.5% natural abundance) to pro- 
duce a potential multiqubit system. We observed Rabi 
oscillations of the 5 = 1/2 and I = Cr spin in a wide 
temperature range of 4 to 290 K and almost indepen- 
dent of field orientation. The measured spin dephasing 
times and ease of material synthesis suggest that such 
spin qubits based on transition metal oxides might be 
suitable for scalable quantum computing. 

Cr doped K 3 NbOs (hereafter Cr:K 3 NbOs) crystallizes 
in a tetragonal unit cell (742m) with lattice parameters a 
= 6.694 A and c = 7.574 A. Fi gure 1 (a) shows a schematic 
of the Cr:K3NbOg system. Nb 5+ ions are shown at the 
corners while Cr 5+ ion is shown at the body center po- 
sition. Single crystals of Cr:K 3 NbOs were prepared as 
described elsewhere [HI E3] and the Cr 5+ concentration 
was determined to be ~ 0.03% in the studied sample. 
Electron paramagnetic resonance (EPR) measurements 
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FIG. 1: (Color online) (a) Crystal structure of Cr 5+ doped 
KaNbOs: Cr [green (or light gray) ball)], Nb [red (or dark 
gray) balls], K (gray balls) and O [blue (or black) balls)]. For 
clarity, the oxygen atoms are shown only around the Cr ion 
(dark gray polyhedron), (b) Room temperature cw EPR spec- 
tra (vertically shifted for clarity) measured at v ~ 9.64 GHz 
for Bo 1 c and -Bo ||c, respectively. In both cases, the strong 
central peak corresponds to the S = 1/2, 1 = resonance and 
the four weak sidebands are due to the 53 Cr isotope (I = 3/2, 
9.5% natural abundance). The inset displays the energy level 
diagram showing the eigenstates of the S = 1/2, 1 = spin 
system in an external field. 



were performed using a Bruker Elexsys 680 spectrometer 
at A-band [y ~ 9.64 GHz) in both continuous-wave (cw) 
and pulsed modes. Temperature was varied between 4 
and 290 K using helium-flow cryostats. 

The electronic spin Hamiltonian for Cr:K3NbOs can 
be written as follows 



H = ix B B Q -g.S + H, 



hf 



H, 



dipole 



(i) 



The first term is the electron Zeeman interaction where 
Hb is the Bohr magneton, g the g-tensor, and S the spin 
operator. Hhf comprises the hyperfine interactions given 
by J2 n S-A-I n . Even though there are no hyperfine inter- 
actions with the 52 Cr (I = 0) nuclei, we have found that 
the superhyperfine interactions from 39 K nuclei (I = 3/2, 
93.3% natural abundance) are non-negligible and also the 
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FIG. 2: (Color online) (a) Hahn echo intensity decay of 
CriKaNbOs as a function of delay time at v ~ 9.64 GHz 
and 25 K for Bo || c. Observed modulation of the echo height 
is due to the superhyperfine coupling with the 39 K nuclei. 
Dashed line is the exponential fit to extract spin-spin relax- 
ation time T%. (b) Temperature dependence of T2 for Bo || c 
at v ~ 9.64 GHz. 



interactions with the 93 Nb (I = 9/2, 100% natural abun- 
dance) cannot be ruled out. The third term represents 
the dipolar interactions between electron spins. In solids, 
dipolar coupling is usually the principle mechanism limit- 
ing the spin-spin relaxation time T 2 . However, in our case 
(~ 0.03% Cr 5+ concentration), a mean separation of ~ 8 
nm between Cr 5+ ions yields an average electron-electron 
dipolar interaction of about 0.1 MHz. Our system is suf- 
ficiently diluted that it can be considered as an almost 
perfectly isolated Kramer's spin system. 

cw EPR spectra of the CnKaNbOg single crystal 
recorded at room temperature are presented in Fig. 1(b). 
The single central peak corresponds to the electron-spin 
transition AS Z — ±1 within the S = 1/2, 1 = doublet, 
as shown in the inset. The g values are calculated as 
g x = 1.9878±0.0002for5 -L c and g { \ = 1.9472±0.0002 
for Bo |j c and are characteristic of a tetragonally dis- 
torted tetrahedral system with a 3d x 2_ y 2 ground state 
(li| . We note that the g anisotropy in our system is small 
and therefore can be treated as quasi-isotropic. The four 
weak satellite peaks flanking the central line for both ori- 
entations arise from the hyperfine coupling to the nuclear 
spin I = 3/2 of the 53 Cr isotope with A\\ = 3.9 ± 0.1 mT 
and A± = 1.1 ± 0.1 mT. In this study, we will address 
only the central resonance from the S = 1/2, / = 52 Cr 
single qubit spin system. 

As a first step, the two important characteristics of a 
qubit namely, the spin-lattice relaxation time (T±) and 
the spin-spin relaxation (T2) time were measured us- 
ing the standard pulse sequences. T± was measured by 
an inversion-recovery [using free induction decay (FID)] 
method employing the sequence 7r-r-7r/2-FID with vary- 
ing t. T\ increases continually from ~ 526 ns at room 
temperature to ~ 1 s at 4 K. This indicates that the spin- 
lattice relaxation is caused by thermal processes. T 2 is 
obtained using a 2-pulse Hahn echo decay sequence 7r/2- 
T-7r-echo that gives the echo intensity as a function of r 
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FIG. 3: (Color online) (a) Cr 5+ free induction decay (FID) at 
4.2 K for two values of applied magnetic field Bo = B r£S + SB, 
with B re3 determined from the cw resonance experiment for 
-Bo || c [see Fig. 1(b)]. Experimental data (squares) are well 
fitted by a damped oscillatory motion (solid lines), (b) The 
beat frequency of the FID versus detuning field. The data 
are well described by a through-origin linear fit. (c) The total 
dephasing time (T2*) as a function of temperature for two 
values of SB [shown in (b) by a triangle and a full circle]. T2 
remains practically constant from 4.2 K to room temperature. 



[see Fig. 2(a)]. The observed oscillations are from the 
electron spin echo envelope modulation (ESEEM) effect. 
Their Fourier transform yields the 39 K nuclear spin levels 
splittings due to the combined effect of nuclear Zeeman, 
superhyperfine and quadrupole interactions. An electron 
nuclear double resonance study at 240 GHz [T3| yields 
the hyperfine couplings which range from 0.41-0.73 MHz 
and which are consistent with the ESEEM results. These 
values allow us to assign the gaussian lineshape and the 
linewidth (4.2 MHz) of the EPR resonance to the unre- 
solved superhyperfine coupling with the surrounding K 
nuclei. For a gaussian line, the relation between the FID 
decay time T 2 an d the linewidth AB pp is given as T 2 
= 2y2/ ABpp. Using the above mentioned AB pp of 4.2 
MHz, T 2 can be estimated as ~ 108 ns. 

The temperature dependence of T 2 for B || c at ~ 9.6 
GHz is displayed in Fig. 2(b). T 2 slowly increases from 
~ 130 ns at room temperature to ~ 10 /its at 70 K and 
essentially remains constant down to 4 K. The details 
will be discussed in a separate paper [ftT ]. 

The decay of the macroscopic magnetic moment per- 
pendicular to the magnetic field due to the dephasing of 
the individual spin packets with respect to each other is 
measured by the FID and is shown in Fig. 3(a). Here, we 
plot the FID measured at 4.2 K for two detuning fields 
5B (= Bq - B res ) for Sn || c-axis orientation. In a rotat- 
ing frame of reference [18|, Bq is along z axis whereas B\ 
is along the x axis. A 7r/2 pulse applied with a +x phase 
rotates the spins, initially aligned along the z axis, to the 
-y axis and when on-resonance (Bq = Bres), the spins re- 
main oriented along the -y axis during the free evolution 
period and decay with T 2 *. When the resonance field 
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FIG. 5: (a) Rabi oscillations damping time tr as a function 
of microwave amplitude at five temperatures for Bo || c ori- 
entation, (b) Rabi oscillations damping rate as a function of 
S}r at 4 K and 100 K for Bo II c orientation. 



FIG. 4: (a) Rabi oscillations are observed as time evolution 
of the (S y ) component as a function of pulse duration at dif- 
ferent power levels and temperatures, (b) Linear dependence 
of the Rabi frequency (Qr) on the microwave amplitude at 4, 
100, 170, 230, and 290 K. (c) Orientation dependence of Q R 
measured at 5 K for two microwave field (Bi) directions. 



is detuned by SB, the average macroscopic magnetiza- 
tion undergoes a circular motion in the transverse plane 
with a beat frequency Sf. The decay time and Sf were 
extracted by fitting the FID data (squares) to a gaus- 
sian damped sinusoidal curve (solid lines). Figure 3(b) 
shows a plot of Sf versus SB. The through-origin linear 
dependence corroborates the relation Sf oc SB. Figure 
3(c) displays the temperature dependence of T 2 * for SB 
— 6 and 8 G. The obtained T 2 * °f ~ H5 ns agrees with 
the value calculated from the observed linewidth and re- 
mains practically constant from room temperature to 4.2 
K. 

We have performed driven coherent spin manipulation 
and obtained the so-called Rabi oscillations. All mea- 
surements were performed by applying a nutation pulse 
around +x axis of length r and measuring the result- 
ing (S y ) by the FID amplitude after a deadtime of ~ 80 
ns. Several examples of the observed Rabi oscillations 
are given in Fig. 4(a) at different power levels and tem- 
peratures. Below 170 K the Rabi amplitude follows the 
expected 1/T dependence. Above 170 K the additional 
decrease of the amplitude (up to 40% at 290 K) is ob- 
served due to the fact that considerable T 2 decay occurs 
during the measurement deadtime. Noticeably, the Rabi 
oscillations are observable even at room temperature. 

The experimental data are well described by a single 
exponential oscillating function 
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The Rabi frequency, £Ir and a damping constant, tr, 
are determined by fitting the data to Eq.(2). In Fig. 
4(b) the Rabi frequency is plotted against the microwave 



amplitude for temperatures of 4, 100, 170, 230, and 
290 K. Since microwave amplitude is proportional to the 
square-root of the incident power, the x axis is taken as 
y/P/Pmax = W A/20 , where P is the power, P max is the 
maximum power (~ 1 kW) and A is the attenuation in 
dB. The linear dependence of Hr on the microwave am- 
plitude is consistent with what is expected for Rabi oscil- 
lations. Since the g value of CriKaNbOs is nearly tem- 
perature independent [lij], the smaller slope of the Qr 
versus incident power at higher temperatures is due to 
the smaller B\ field at the sample caused by the smaller 
resonator Q-factor at higher temperatures. The depen- 
dence of Qr on the direction of the microwave field, B±, 
at 5 K is displayed in Fig. 4(c). We note that ratio of 
£Ir between B\ || c and B\ _L c, (CIr n) /(Oh,x)) deviates 
from the expected g\\/g± ratio and is likely due to the 
small changes in the resonator Q factor with the sample 
orientation. However, this small directional dependence 
of £Ir on the microwave field confirms the quasi-isotropic 
nature of our system. In addition, we have also checked 
the directional dependence of FIr on the external field, 
Bq and find no difference between the two orientations 
(not shown here). 

The variation of Rabi decay time (tr) with drive am- 
plitude at different temperatures is shown in Fig. 5(a). 
With decreasing microwave power tr increases and then 
saturates in the low power limit regardless of tempera- 
ture. In addition, the saturation value of tr increases 
strongly upon cooling from 290 K to 100 K and then 
remains more or less constant between 100 K and 4 K. 
This behavior resembles the temperature dependence of 
T 2 [see Fig. 2(b)], indicating that tr is restricted by a 



T 2 mechanism. As in ref. 14j this can be explained by 
a model based on a stochastic field with a normalized 
distribution (3 [HHo)], 



1/tr = f3n R + 1/2T 2 



(3) 



In our case, this stochastic field can arise from the mod- 
ulation of the superhyperfme couplings and the radiation 
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induced changes in the dipolar field [20| . We find that at 
high temperatures tr ~ ITi as expected. In the low tem- 
perature range a marked dependence of the Rabi oscilla- 
tions decay rate on the Rabi frequency or the amplitude 
of the B\ field is observed. However, this dependence 
[see Fig. 5(b)] does not show the simple linear depen- 
dence from Eq.(3) which was observed in AIO4 - and E' 



centers in silicate |20|, l21| , and extrapolates to about 2 fj,s 
at temperatures below 100 K, which is somewhat shorter 
than the observed T2 (2.7 — 10 /is) in this range [see Fig. 
2(b)]. We tentatively ascribe this to the larger role played 
by the 39 K superhyperfine interactions in our system. 

The single qubit figure of merit Qm , defined as Qm = 
QrT2/tt [22j, estimates the efficiency of a quantum de- 
vice. Using a T2 time of 10 /is we obtain ~ 500 coherent 
single-qubit operations at liquid helium temperature, im- 
plying that Cr 5+ doped KaNbOs is a viable electron spin 
qubit. One way to improve the Qm is to increase the B\ 
by going to higher microwave power. Even though the T2 
of Cr:K3NbC>8 is smaller than that of nitrogen-vacancy 
centers in diamond [9L [icil| . it is comparable to that of 
the rare-earth qubits [l4 | and molecular magnets [l3j |. 
As mentioned earlier, the spin decoherence in our system 
is likely due to the superhyperfine interactions with 39 K 
nuclei. This hints a possibility of synthesizing samples 
with even longer coherence time. For example, replace- 
ment of 39 K by 41 K reduces the superhyperfine coupling 
by nearly a factor of 2. We have also succeeded in synthe- 
sizing KCaCrOs, which should reduce the superhyperfine 
interactions due to 39 K by two thirds. 

Qubits could be coupled either by dipolar (as in Q) 
or by exchange interactions that are of the order of few 
Kelvin in this system 23j . We have observed that at con- 
centrations higher than 5%, we do form Cr-Cr exchange- 
coupled triplet pairs, which should in principle enable 
clustering of several qubits. Such clusters could be fur- 
ther coupled either by dipolar or exchange interactions. 

For the present studies, sub-mm sized single crys- 
tals were synthesized from solution. In future studies, 
this technique will be modified to allow precipitation of 
nanocrystals in areas designed by lithography on an elec- 
tronic circuit. The qubits could be entangled via photons 
when integrated in microcavities or via current oscilla- 
tions if spin detection is done by electronic transport. It 
is important to note that, in contrast to the 2DEG-based 
quantum dots where the free electron overlaps over the 
dot's nuclei, the Cr spins interact only with the neigh- 
boring non-zero spin nuclei of the crystal. Spin detection 
can be done in several ways, for instance by optical de- 
tection 24 1 , single-electron transistors [25| , nanoSQUIDs 
[2^ ] , or microwave detection by means of SQUID [27| or 
Hall probes [3]. The quasi- isotropic magnetic character 
of the Cr spins demonstrated here is relevant for most 
practical implementations when a controlled positioning 



of the spins is hard to achieve. Therefore, we consider 
this novel material to be potentially well suited for inte- 
gration with standard nanofabrication methods used for 
on-chip studies. 

In conclusion, we have reported on the observation of 
Rabi oscillations over 4 — 290 K in CriKaNbOs, an es- 
sentially pure S = 1/2, metal-oxide system. We find the 
intrinsic phase-coherence time T2 ~ 10 ^s and the single 
qubit figure of merit Qm ~ 500 at liquid helium tem- 
perature. Our results demonstrate that the transition 
metal oxide-based spin systems hold high potential for 
quantum information applications. 
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